Textures and the shapes of domains in Langmuir monolayers 
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Two-dimensional domains containing an XK-like order parameter exhibit non-trivial internal struc- 
ture and take on shapes controlled by the configuration that the order parameter adopts. The 
textures exhibited by the order parameter in such domains are controlled by the interplay between 
bulk and surface contributions to the energy. We report calculations of the internal texture and the 
shape of such domains. These calculations lead to the determination of the equilibrium properties 
of two-dimensional domains, such as those observed in Langmuir monolayers. This allows for the 
unambiguous exploration of the implications of experimental findings. 
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Monolayers of surfactants confined to the air/water in- 
terface have been found to possess complex textures sim- 
ilar to those observed in liquid crystals. The textures are 
generally observed in "tilted" phases — that is, in phases 
in which the long axes of the molecules in the film are 
not perpendicular to the water surface but are uniformly 
tilted with respect to the normal. The textures are the 
result of the spontaneous organization of the molecular 
tilt azimuth on macroscopic length scales. They can be 
understood, at least qualitatively, in terms of continuum 
elastic theories of smectic liquid crystals or, more rele- 
vant to the present work, in terms of the two-dimensional 
XY model. 

When the surfactants organize into domains of con- 
densed tilted phases, such as the Li phase, surrounded by 
an isotropic phase 0] , many striking textures have been 
observed. Such textures range from continuous boojums, 
textures with point and/or line defects Q, to striped 
textures in spiral domains We will, however, focus 
our discussions on the boojums in which the tilt azimuth 
varies continuously and appears to radiate in some cases 
from a defect located at the edge of the domain or from 
a "virtual" defect in the isotropic phase 0]. Domains 
containing a boojum texture are not circular. Rather, 
they reflect the various nonisotropic influences that con- 
trol domain shape in the presence of a nonhomogeneous 
background. Among the features seen are protrusions, 
at times sharp enough to be characterized as "cusps" [0 , 
and indentations, which give the domains a heartlike, or 
"cardioid" appearance (|J^]. 

While there are notable cases in which the full analyti- 
cal solution can be obtained , the determination of the 
properties of a domain in which such a texture has formed 
presents a calculational challenge. Energy minimization 
entails the simultaneous adjustment of the orientations 
of surfactant molecules in the interior of the domain and 
the shape of the domain's boundary. Past work on the 
problem of the shape of surfactant domains has relied on 
either the assumption of a texture that is unaffected by 
variations in the boundary Isj, or on the approximation 



of a nearly circular domain |3||9|. Neither assumption is 
necessarily close to reality, and results obtained in both 
cases can be legitimately called into question. 

This Letter describes a successful evaluation of equi- 
librium properties of a domain of surfactants. The evalu- 
ation is based on the numerical solution of the extremum 
equations for the energy of an XY-like texture confined 
to a compact, but not necessarily circular, domain. We 
are able to explore the texture and the associated domain 
shape that result from various forms of the boundary and 
bulk energy of the two-dimensional system. 

The starting point in the analysis is the Hamiltonian 
of the XY order parameter, represented as the two- 
dimensional unit vector c(x,y), which is parameterized 
in terms of the angle, Q(x,y) between c(x,y) and the 
x axis. The order parameter field @(x,y) will be re- 
ferred to as the texture. The Hamiltonian, H[Q] — 
f n HbdA + <f r <r{d — 9) ds, where Tib is given by 



H b = U\VQ\ 
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<d x and <d y represent partial derivatives of 9 with re- 
spect to x and y, respectively. The integral J n is over 
the domain's bulk, while j> r is an integral over the closed 
curve bounding the domain. The coefficient k is the 
"mean" Frank constant, the average of the bend and 
splay moduli, while (3 is proportional to the difference 
between the two moduli. Specifically, 2k — K s + Kb 
and 2kP = K s — Kb where Kb/2 multiplies |V x c| 2 in 
the energy of the XY order parameter and K s /2 mul- 
tiplies |V • c| 2 . The constants Kb and K s are, respec- 
tively, the bend and splay modulus. The boundary en- 
ergy, a ($ — 9) , will have the general form 



cr(ip) = ao + a\ cos ip + et2 cos 2(p + • 



(2) 



The angle & in the argument of the boundary energy is 
the angle between the unit normal to the curve bounding 
the domain and the x axis. The fact that the harmonic 
expansion, (H), of the boundary energy consists entirely 
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of cosine terms reflects the absence of "chiral" interac- 
tions between the texture and the domain boundary. 

The minimization of the energy leads to equations for 
the texture <d(x,y) and the bounding curve T. 6(2;, y) 
satisfies both a bulk and a boundary extremum equation. 
The bulk equation is 

- V 2 6 + /3 [(Q xx - Q yv ) cos 26 + 2Q xy sin 26 

+ + Q 2 y ) sin 26 + 26^6^ cos 26] = 0. (3) 

The double-subscripted 6's represent second partial 
derivatives with respect to the relevant variables. The 
extremum equation for Q(x, y) on the boundary is 

K6„[l-/3cos2(?9-e)] 
+K/36 t sin 2(^-6) - </(i? - 0)= 0, (4) 

where 6„ and 6t are respectively the normal and tangen- 
tial derivatives. The extremum equation for the bound- 
ing curve r is 

n b - - 6)e n - o"{4 - e)e t 

d-d 

+ [a(#-Q) + a"(ti-e)} — + A=0. (5) 

The quantity A is a Lagrange multiplier that enforces the 
condition of constant enclosed area. 

We first solve numerically for the equilibrium Q(x,y) 
using a variational formulation of the finite element 
method |l(|, assuming an initial boundary . A mesh 
of triangles is generated over f2 using an adaptive method 
that refines the grids where necessary. Functions are de- 
fined by their values on the vertices of the triangles. The 
value of a function elsewhere is obtained by interpola- 
tion. The system energy H [6] is then a function of the 
values of Q(x, y) at the vertices, = (6A which can be 
determined by solving dH(Q)/d®i = fill ; where i runs 
from 1 to the number of vertices. We label the texture 
determined at this point 0(°). 

Equation (H), after being cast into a coordinate- 
dependent form, turns out to be a second order differen- 
tial equation for T. Substituting ®(°> for the Q(x,y) in 
the equation, T can be determined by the Runge Kutta 
method JTlj] . The solution, which is labeled is in 

turn utilized to determine a new texture The pro- 

cess is iterated until self-consistency is achieved. 

Making use of the numerical scheme outlined above, 
we have investigated the the texture and shape of the 
bounding curve under the influence of variations in the 
boundary energy coefficient a,2 in Eq. (||) and the stiffness 
coefficient (3. Before presenting our results, we note that 
when et2 = and (3 — 0, the exact result Q is given 
by a circular boundary of radius R, with the boojum 
texture, i.e., a defect with winding number +2 || located 
at a distance Rb = R(l + \/l + P 2 )/p from the center of 
the domain, where p = Ra\/ k is the normalized domain 
radius. The normalized domain radius p will be used as 




FIG. 1. The constant-order-parameter contours and shapes 
of domains computed for —0.5 < 02 < 0.5, p = 5, do = 4, 
ai = 1. (a)oa = -0.5, (b)o a = -0.3, (c) a 2 = -0.1, ( d) 
a-2 — 0.1, (e) a2 = 0.3 and (f) 02 = 0.5. 




FIG. 2. The shapes of domains of various sizes computed 
for ao = 4, a\ = 1 and 0,2 = 0.6. (a)Smaller domains 
with p = 0.2, 0.25, 0.33, 0.5, 1 which exhibit 2-fold symme- 
try. (b)Larger domains with p — 2, 2.5, 3.3, 5, 10 which have 
a protrusion on one end of the boundary. Each of the sets of 
domains are plotted to scale. 

a gauge of the domain size throughout the discussion. 
For cases in which the domains are not circular, R is 
the effective radius, in that ttR 2 will be the area of the 
domain. 

We first look at our result by varying only 02 while 
keeping (3 = 0. In the limit of very small domains, p <C 1, 
the line tension anisotropy has very little effect on the 
texture, which exhibits almost no spatial variation. The 
boundary response of the domain is not affected by the 
a\ contribution. If 02 ^ 0, small domains exhibit twofold 
symmetry. They become elongated when a<i > 0, and 
they flatten at both ends when 02 < 0. For larger do- 
mains in which p ~ 1 and 02/01 <C 1, we find a pro- 
trusion when d2 > 0, and an indentation when 02 < 0. 
The feature lies on the end of the domain boundary clos- 
est to the virtual defect. The loss of twofold symmetry 
is largely due to the fact the dominant contribution of 
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FIG. 3. The constant-order-parameter contours and the 
shapes of domains with p = 1 and ao = 4. The coeffi- 
cients of the anisotropic line-tension are (a)ai = 1, 02 = 0, 
(b)ai = 0.8, a 2 = 0.2,(c)oi = 0.6, a 2 = 0.4,(d)ai = 0.4, 
(J2 = 0.6,(e)oi = 0.2, (J2 = 0.8 and (f)ai = 0, 0,2 = 1. 

the domain texture, which in turn induces deformation 
through a 2 , comes from the a± contribution in the bound- 
ary anisotropy. Boundary features are more prominent as 
the domain size increases. Figure |l| shows the shapes of 
domains for —0.5 < a 2 < 0.5 for p = 5. Numerically, the 
boundary features are well-defined for domains of sizes 
p ~ 1. To the numerical accuracy that we are able to 
achieve, there is no discontinuity in slope on the domain 
boundary. That this ought to be the case can be veri- 
fied analytically |IJ,[l3|]. The numerical result we have 
is compatible with those discussed in Refs. For 
even larger domains, p 1, the features are confined 
in a small portion of the boundary. The domains be- 
come nearly circular again in the large-/? limit. Up to the 
largest domain we have examined, p — 32, we can iden- 
tify protrusions when 02/00 > 0.1. Because of the rapid 
texture variation in the immediate vicinity of the bound- 
ary, associated with the approach to the boundary of the 
virtual boojum singularity, we are unable to perform de- 
pendable numerical investigations of extremely large do- 
mains. This leaves open the question of the asymptotic 
behavior of the domain in the large-/? regime. Figure ^ 
shows the shapes of domains in which the normalized 
radius ranges from p = 0.2 to p = 10 for a\ — 1 and 
a 2 = 0.6. 

As regards the texture, we are able to numerically re- 
produce the exact result of Rudnick and Bruinsma ||] 
when eti = and a 2 = 1 and the boundary is fixed at a 
circle. This texture is associated with two virtual defects. 
When the boundary is allowed to relax, the domain ac- 
quires a "cigar shape" . When both a\ and a 2 are not 



FIG. 4. The constant-order-parameter contours and the 
shapes of domains with p — 8, ao — 4, a\ — 1.6 and 02 = 0. 
Their stiffness coefficients are (a)/3 = —0.8, (b)/3 = 0.8. 




FIG. 5. The domain shapes computed for ao = 4, aj. = 1.6 
and 02 = 0, and p — 0.5, 1, 2, 4, 8. Their stiffness coefficients 
are (a)/3 = —0.8, (b)/3 = 0.8. For ease of observation, domains 
are not shown to scale. 

equal to zero, the texture can be thought of as a super- 
position of pure a\ and pure a 2 textures. Typically for 
p ~ 1, the effect of the second defect becomes observable, 
in the form of a distortion of constant-order-paramcter 
contours, when a 2 ja\ w 1/4. The progressive changes of 
the texture and domain shape from a system with pure 
a\ to a system with pure a 2 are depicted in Fig. |^. Do- 
mains with indentations, protrusions, and cigar-shaped 
domains, have all been observed (j. 

We now investigate the textures and boundaries that 
result when 0^0. The parameters ao and a\ are as- 
sumed to be finite, while all other a$'s are set equal to 
zero. We have computed the textures and domain shapes 
for -0.8 < < 0.8 and 0.5 < p < 8. When \(3\ < 0.5, 
our results are in qualitative agreement with those re- 
ported in Ref. ||. For < 0, the texture is modified 
as if the virtual singularity has moved closer to the do- 
main boundary, and the domain acquires a very small 
protruding feature. On the other hand, when > 0, the 
texture relaxes as if the virtual defect has moved away 
from the boundary, and the boundary correction is simi- 
lar to that of a domain with an indentation. There are, 
however, quantitative differences. The textures obtained 
numerically deviate significantly from the boojum tex- 
tures utilized in Refs. |@. The sizes of the features on 
the boundaries are no more than 1% of the overall domain 
radii. We do not expect such features to be observable 
experimentally. 

When |/3 1 = 0.8 and p = 8, the domains are no longer 
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circular in appearance. Domains acquire small protru- 
sions for negative (3. Domains have an indentation when 
[3 is positive. Figure || shows the domain shapes and tex- 
tures computed with (3 — ±0.8 for p = 8. The deviation 
from the boojum texture is more clearly observable when 
[3 = +0.8 in Fig. |]b. The constant-order-parameter con- 
tours distort so as to push the virtual defect away from 
the boundary. Figure || shows the changes in the domain 
shape when p decreases from 8 to 0.5 for [3 = ±0.8. The 
amplitude of the protrusion of the domain with nega- 
tive [3 decreases with p. At the same time, segments of 
the boundary with negative curvature start to develop 
adjacent to the protrusion. At p = 2, the negative curva- 
ture segments of the boundary overwhelm the protrusion 
and give the domain an indentation. When p = 0.5, the 
indentation disappears, while the protrusion continues 
to decrease in amplitude. The domain appears circular. 
The complex boundary behavior described above can be 
attributed to the rapid variation of the order parameter 
along the boundary near the singularity. In contrast to 
this behavior, the size of the indentation on the bound- 
ary decreases monotonically as p decreases when (3 = 0.8. 
The domain appears circular at p < 2. It is the strength 
of the splay modulus that smooths out the variation of 
the order parameter along the segment of the boundary 
near the singularity and results in simple boundary be- 
havior. 

Our attempts to compare our results with experiments 
indicate that all the observed domain shapes in Refs. 
|7|,p|JT4]| can be accounted for, a least qualitatively, by 
line-tension anisotropy alone. Although we have estab- 
lished the existence of nontrivial domain shapes due to 
purely elastic anisotropy, domains with the protrusions 
and indentations discussed immediately above occur only 
when the anisotropy is very strong \f3\ — 0.8. Fur- 
thermore, the protrusions that are generated in our cal- 
culations do not resemble those reported to have been 
observed experimentally pjl^l- We conclude that line- 
tension anisotropy must be present in such monolay- 
ers. On the other hand, both line-tension and elastic 
anisotropics may be responsible for indentations on the 
boundary of domains of the compressed phase. Fur- 
ther measurements on the size dependence of the domain 
shapes are needed to identify the underlying mechanism. 

We note here that our analysis is based on a highly 
simplified model of monolayers. Other factors that con- 
tribute to the boundary shapes and textures — such as 
the dipolar interaction the tilt degree of freedom 
[M and the bond-orientational ordering [Q — have been 
neglected. However, the results reported here are useful 
in that they do indicate the influence of various features 
that make up the model discussed. 

In conclusion, we have devised a numerical scheme that 
enables us to solve simultaneously for the minimum en- 
ergy configuration of an XY order parameter confined to 
a two-dimensional domain and the extremal shape of the 



boundary of that domain. We have utilized the method 
to investigate the response of the texture-boundary sys- 
tem under the influence of a nontrivial anisotropy in the 
boundary energy as well as the bulk elastic energy. Both 
line-tension anisotropy and elastic anisotropy result in 
nontrivial domain shapes. The observed domains with 
protrusions and cigar shape can be attributed to the line- 
tension anisotropy. On the other hand, elastic anisotropy 
cannot be ruled out as the source of the observed car- 
dioids. 

We are grateful to Professor Charles Knobler, Profes- 
sor Robijn Bruinsma and Dr. Jiyu Fang for very useful 
discussions. 
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